ABSTRACT 



We investigate the magnetic field behaviour of an antiferromagnetic Heisenberg 
spin-1 chain with the most general single-ion anisotropy. We discuss the regime in 
which the magnetic field is below the transition value. The splitting of the Haldane 
triplet is obtained as a function of a field applied in an arbitrary orientation by 
means of a Lanczos exact diagonalization of chains of up to 16 spins. Our results are 
nicely summarized in terms of a first-order perturbation theory. We explain various 
level crossings that occur by the existence of discrete symmetries. A discussion 
is given of the electron spin resonance and neutron scattering experiments on the 
compound Ni(C2H8N2)2N02C104 (NENP). 
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I. INTRODUCTION 

There is good evidence that the integer-spin Heisenberg antiferromagnetic chains have 
a gap, as suggested by Haldane^. This is estabhshed on the theoretical side by various 
techniques including finite-size calculations as well as field theoretic arguments^. On the 
experimental side the first evidence^ came from CsNiCla: neutron scattering (NS) re- 
vealed an excitation gap. Since this compound is only moderately one-dimensional, para- 
sitic three-dimensional effects complicate the picture. On the other hand, the compound^ 
Ni(C2H8N2)2N02C104, abbreviated NENP, is much more one-dimensional and remains in 
a magnetically disordered state even at very low temperature. Zero-field neutron scattering 
experiments^ have clearly shown the existence of the Haldane triplet split by easy-plane 
anisotropy. In addition, a number of other measurements have been performed^"^ in a 
magnetic field: susceptibility, high-field magnetization, and neutron scattering in a finite 
field. The application of a magnetic field leads to a Zeeman splitting of the Haldane triplet 
and one member of this triplet crosses the ground state at a critical value Hg that depends 
on the field orientation. This is clearly seen in NS experiments where all members of the 
triplet can be followed individually^. Experiments using electron spin resonance technique 
(ESR) are in excellent agreement with NS as is the case for far-infrared spectroscopy 
measurements ■^■'^ and there is at the present time a satisfactory picture of the behaviour of 
NENP from the experimental point of view. 

On the theoretical side, effective quantum field theories have been used to predict 
the magnetic field behaviour of the spin-1 chain Heisenberg Hamiltonian that models the 
magnetic properties of NENP. In fact the original work of Haldane showed that, in the 
large integer spin limit, the antiferromagnetic spin chain is described in the low-energy 
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limit by an 0(3) nonlinear sigma model. This nontrivial field theory is difficult to study 
and, based on the large-N limit of the 0(N) model, it has been suggested that a simple 
theory with three massive bosonic fields might be an appropriate effective theory of the 
spin-1 chain. It is then possible to obtain the behaviour^^'^^ of the system in a magnetic 
field. 

Another possibility has been suggested starting from an integrable chain It is known 
that the spin-1 Hamiltonian Hq — Sj • Sj_|_i — (Sj • Sj_|_i)^ is solvable by the Bethe 
Ansatz technique^^ and leads to a massless theory. This massless theory is characterized 
by an SU(2);i;=2 symmetry and can be realized by three (massless) Majorana fermions. 
As one perturbs the Hamiltonian Hq towards the pure Heisenberg Hamiltonian without 
biquadratic coupling, it is natural in the framework of the Haldane conjecture to expect 
that these fermions become massive. One can then use a theory of three free massive 
Majorana fermions to approximate the Heisenberg chain. These two theoretical approaches 
are not in complete agreement and it is thus interesting to have results of a completely 
different nature, based on numerical studies of finite chain diagonalization. 

In this paper, we present the results of our study of the field behaviour of a spin-1 
chain including realistic single-ion anisotropies of the most general kind. We diagonalize 
by means of a Lanczos algorithm chains of up to 16 spins under a magnetic field applied 
in various positions. Our findings are neatly summarized by a simple perturbation calcu- 
lation that may be used as a practical tool to obtain the field behaviour of an Haldane 
magnet. When the applied field become strong enough there is a phase transition towards 
a magnetically ordered phase^^'^^. In this paper we will restrict ourselves to the singlet 
phase where the Haldane gap is not destroyed. We give a discussion of the various level 
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crossings that may or may not appear depending on the field orientation with respect to 
the symmetry axis of the crystal. Section II contains the treatment of the magnetic field 
as a perturbation. Section III explains the Lanczos results as well as their relationships 
with the perturbative expansion. Section IV contains our conclusions and a discussion of 
NENP experiments. 

II. THE PERTURBATIVE RESULTS 

We focus on the microscopic Hamiltonian of a spin-1 antiferromagnetic chain in an 
applied field and single-ion anisotropy: 

7i = J ^ S, . S,+i + D{Sff + E[{Sff - {Sff] -n Si. {II - 1) 

i 

Here Sj are quantum spins S=l and we include the Bohr magneton and the Lande g factors 
in the definition of the magnetic field. The exchange coupling J is taken to be positive i.e. 
antiferromagnetic. Below we set J — 1. The D and E terms in Eq.(II-l) parametrize the 
most general single- ion anisotropy. We work with periodic boundary conditions. In zero 
field and in the absence of anisotropy [D = E = Q) the Hamiltonian (II- 1) is invariant 
under the full SU(2) rotation group. If the D term is nonzero then the symmetry is broken 
to a residual U(l) subgroup of rotations around the Z axis. If, in addition to D, there is 
some further in- plane anisotropy (nonzero E) even this U(l) is broken. However in this 
case there are still discrete remnants of the initial SU(2): the system is invariant under tt 
rotations around the coordinate axis X, Y, Z. These symmetry operations are denoted by 
RZ,Ry, Rz in this article. 

If we now add a magnetic field, there is further symmetry reduction: when H is in a 
generic position with respect to the coordinate axis, even the discrete operations R'^,Ry, R^ 
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are lost. However, there are some special orientations of H that retain discrete symmetries: 
if H lies along the a axis the symmetry is conserved. This corresponds simply to 
rotation around the magnetic field axis. But there are also discrete symmetries when the 
magnetic field lies in one of the planes {X, Y), (F, Z) or {X, Z). This is seen easily when 
the Y-component of H is zero: the Hamiltonian expressed in the basis of eigenstates of Sf 
is then a real symmetric matrix. Complex conjugation is thus a discrete symmetry. In fact 
Complex conjugation simply changes the sign of the operators Sf and is a symmetry when 
Hy = (in this particular basis). When the magnetic field is in another plane the system 
is still invariant under a combined operation involving complex conjugation and a rotation 
K^. If, for example, Hg = then complex conjugation will transform HySf into —HySf 
and then one uses i?J to come back to the original Hamiltonian. These extra symmetries 
will explain the various level crossings that are found in our exact diagonalization studies 
reported in section III. 

We now discuss the application of perturbation theory with respect to the D, E, and 
H terms in Hamiltonian (II-l). We write Eq.(II-l) as: 



The Hamiltonian Tio has full rotational symmetry and its levels can thus be classified 
according to their spin. We focus on the effect of the perturbation on a singlet and a triplet 
state. We know that the ground state of TIq is in fact a singlet with chain momentum 
K = and the first excited state is a triplet with momentum K = tt. Our statements about 
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TL = "Ho + Ti,!, 




{II - 2) 



the effect of the perturbation are general since they are dictated by the Wigner-Eckhardt 
theorem. 

Let us first discuss the case of a singlet state |0). Due to complete isotropy, the 
following equalities hold: 

(0| J2iStr\0) = (0| ^(-S]^)^|0) = (0| J2{Stf\0) = InS{S + 1) = ^N. {II - 3) 

i i i 

In addition the vector (0| Ylii Si|0) is zero. We thus obtain the first order shift of the singlet 
energy: 

= Eo + D{Q\ Y,{Stf\^) ^Eo + ^DN. {II - 4) 

i 

We note that there is no eff'ect coming from the in-plane E-term or H. We now discuss 
the triplet splitting in zero field. In the standard basis the triplet states are noted |lm), 
m = —1, 0, +1. Perturbation theory involves matrix elements of the following operator: 



Qaf3 ^ J2 SfS^ - {II - 5) 



It is a spin-2 irreducible tensor operator since it transforms as a traceless symmetric ten- 
sor. Its standard components are noted M = -2, . . . , +2. The Wigner-Eckhardt 
theorem implies that the matrix elements in a triplet state of such an operator are related 
by: 

(lm|C»(2^)|W) = C ■ (lm|21Mm'). (// - 6) 

In this equation (lm|21Mm') is the Clebsch-Gordan coefficient coupling two spin-1 states 
to a spin-2 state. Thus perturbation theory is characterized by a single coefficient C (to 
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first order). This is easily found when working out matrix elements in the canonical basis 
\x), \y), \z) rather than using the standard basis. We consider the matrix elements: 

i 

Application of a tt rotation shows that this element is zero if a and 7 are distinct. In this 
case there are only two different matrix elements. We can choose: 

{x\J2{S!r\x){^a) and {x\J2{Sff\x){^b). {II -8) 

i i 

In addition one notes that b+2a = NS{S+1) = 2N. The perturbation Hs = Ei + 
-^[("^f )^ ~ i'^i)'^] ^he canonical basis for the triplet can thus be written as: 

/aD + {2N -3a)E ° \ 

ns= I aD + {3a- 2N)E . (// - 9) 

\ 2{N-a)Dj 

As dictated by Wigner-Eckhardt theorem there is only one coefficient that characterizes 

the perturbation. Subtracting the ground state energy one is led to the following gap 

values: 

= A - kD + 3kE 
Ay = A-KD- 3kE {II - 10) 

= A + 2kD. 

In this equation A is the unperturbed triplet-singlet gap and we note k, = 2N/3 — a. The 
coefficient a is extensive but the difference 2N/3 — a is finite in the thermodynamic limit 
since k appears in gap values. The value of k is not dictated by rotational symmetry of 
course and its numerical value depends for example of the moment of the triplet. When 
E — Q this splitting has been studied in detail^^'^^ for the Haldane triplet with chain 
momentum tt which is the lowest lying triplet. We note that already for D ^ 0.2J there 
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are deviations from the previous perturbation theory: the slopes with D of the two gaps 
are found to be: 

= Aj/ = A - 0.57i:> 

(11-11) 

A^ = A + 1.41i:>. 

This is the best hnear fit of the Lanczos results between D = 0.10 J and D = 0.25 J. The 

slope ratio is already slightly different from 2 as given by Eq. (11-10). The curvature of the 

gaps as functions of D is clearly seen in the data of ref.[17]. In fact a fit including quadratic 

terms in D of the same data leads to the following result: 

Aa; = Ay = A - 0.668i:' + 0.269L>2 

(// - 12) 

A^ = A + 1.357D + O.UbD^. 
The first-order terms satisfy the perturbative results (11-10) and we clearly see the deviation 

from first-order. For practical purposs it is simpler to use the fits (11-11). 

We now add the magnetic field to the perturbative treatment. The ground state, 

being a singlet, is not affected at first order and only the triplet changes. The matrix of 

the perturbation in the canonical basis is then: 

(Px iHz -iHy \ 
-iH, Py iH^ . (// - 13) 

iHy -iH^ p, J 

The quantities Pa denote the diagonal matrix elements of Eq.(II-9). Shifting the origin 
of the energies it is convenient to set pa = Aq, and then the eigenenergies of (11-13) are 
directly the gaps. In section III we show that the resulting values are always extremely close 
of the Lanczos results when the magnetic field is below the critical value. Diagonalization 
of (11-13) does not lead to compact formulas except when the field lies along one of the 
symmetry axis. If H lies along the direction 7 we denote its only nonzero component by 
H^. The eigenvalue A^ is unperturbed and the two other eigenvalues are given by: 



1 
2 



(// - 14) 



Here a and /? are the two other coordinates. There is hyperbohc repulsion of the two gaps 
A+ and A~. The smallest gap A~ goes to zero for a critical value = Aq,A^. At large 
fields the asymptotic behaviour of the gaps is linear. Of course the crossing of one member 
of the triplet with the ground state signals a phase transition^^'^^ beyond which we do not 
expect to gain something from a simple perturbative calculation since other states with 
higher spins also cross the ground state above Re- 

It is interesting to note that this hyperbolic behaviour (11-14) is exactly what is found 
in the fermionic effective theory-^^ for states with K = 'k. When the field becomes large 
and parallel to a coordinate axis the eigenstates of (11-13) take a simple form: they are 
given by la) ± i\fi). We note that the vanishing of A~ at He occurs linearly contrary to 
the free boson prediction. 

If the gap value A^ (which does not move with the field) lies above or below the 
two gaps Aq, and A^ there is in general a crossing of levels between one of the gaps A='= 
and A-y before the critical field. As shown in sect. Ill since they behave differently under 
exact discrete symmetries of the system, we expect that these crossings will survive beyond 
perturbation theory. If however the magnetic field is no longer in a high symmetry position 
these symmetries are broken and one should see only avoided crossings. 

III. LANCZOS RESULTS 

We have performed a Lanczos study of the Hamiltonian (II- 1) on chains of lengths 
N=4,6,8,10,12,14,16. For a generic orientation of the magnetic field, there are no symme- 
tries available apart translational symmetry to reduce the size of the problem. Thus many 
iterations were required to get the first few excited levels. The energies of the ground state 
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(in the subspace K = 0) and the three low-lying levels (in the subspace K = tt) have been 
obtained with a typical precision of 10~^. For a generic orientation of the magnetic field, 
the size of the complex Hermitian matrix H, 3^ , is reduced by translational symmetry to 
fa 3^ /N. The size of the Hilbert space for N=16 is ~ 1.3 x 10^. On one Cray-2 processor, 
acting with H on a vector takes about 17 seconds and the precision of 10~^ is reached for 
fs 60 iterations. 

We have followed the triplet splitting as a function of the E-term. Neutron scattering 
experiments-*^^ have shown that there is a small splitting of the two low-lying modes in 
the case of NENP: A^; « 1.05 meV and 1.25 meV. This means that the E-term is 

much smaller than the D-term in Eq.(II-l) 2.5mey. We treat it perturbatively 

but keep the D-term in the unperturbed Hamiltonian to be solved by Lanczos technique. 
When E=0 the Haldane triplet is split in a high-energy singlet and a doublet. First-order 
perturbation theory for the E-term requires its matrix elements in the subspace spanned 
by the doublet. We find a linear splitting: 

{D) + Ko{D)E 

(III - 1) 

Ay = Ao{D) - Ko{D)E 
Here Ao(-D) is the doublet gap when E=0. It is known by eq.(II-ll). The constant ko{D) 
has been computed for all lattice sizes. Its values as a function of D are given in fig.l. We 
have used the Shanks algorithm^^ to obtain an estimate of the thermodynamic limit value 
of K. In the case of NENP D Ri 0.18 we find in the thermodynamic limit kq ~ 2. The 
corresponding estimate for the in- plane anisotropy is £^ 0.012. 

We now take for granted the zero-field splittings and add a magnetic field along the 
coordinate axis. We present the Lanczos results for the longest chain we were able to deal 

10 



with in figs. 2-6. We find that the numerical points are very well reproduced for all chain 
lengths by the following procedure: we take the zero-field gaps as inputs in the perturbative 
formula (11-14) and obtain the field behaviour. The corresponding curves are plotted as 
solid lines in figs. 2-6. One has of course to vary the gaps with the chain length but the 
hyperbolic behaviour holds for all lengths. The deviations between the exact results and 
the perturbative curves are of the order of 10~^. In figs 2,3,4 the field lies respectively 
along Z, X and Y. One of the gaps is barely affected and the other two are split according 
to the simple formula (11-14). If H lies along Z there is a crossing between the and Aa; 
modes (fig. 2), if H lies along Y then there is a crossing between A^; and A^^ (fig-4). 

It is interesting to note that, before the critical field is reached, one observes that 
other states with higher spin begin to arrive from higher energies. This is seen in all our 
figures: the upper "triplet" mode always follow the perturbative trend but there are new 
states that are lower in energy close to He that do not belong to the Haldane triplet. These 
states will ultimately cross the ground state after the transition point^^. 

Let us consider the case of fig. 2 where H lies along Z. Then the operation i?J is a 
symmetry operation of the Hamiltonian. We can classify the triplet members according to 
their behaviour under K^: 

r:\z)^+\z), {HI -I) 

while 

Rl\x) = -\x) and Rl\y) = -\y). {HI -2) 

Thus no matrix element can avoid the crossing between these two members of the triplet 
(x and z). It is important to realize the following: when H = and Tis = we can 
demonstrate (111-1,2) for the degenerate triplet. By continuity the behaviour under the 
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preserved symmetry operations will survive the addition of H and Hs- This does not 
rely upon perturbation theory but rather on the continuity of a discrete quantum number 
as a function of field and anisotropy. This reasoning holds also for X and Y axis. 

When the field no longer lies along a coordinate axis, the crossings are avoided. A 
typical example is given in fig. 5. Here the field is very close to the Z axis and thus 
there is an avoided crossing between the upper mode (z) and the intermediate (x) mode. 
This should be compared with fig. 2 (H//Z). The solid lines in fig. 5 are obtained from 
the diagonalization of perturbation (11-13). It is always very close to our Lanczos results. 
We have checked that for various field arrangements the same property is true. Another 
example is given in fig. 6 where the field lies close to the Y axis. 

IV. CONCLUSION 

We have studied the magnetic field behaviour of a realistic spin-1 chain with the most 
general single-ion anisotropy that has proved adequate for NENP. The three gaps have 
been computed as a function of the field by means of a Lanczos technique for chains of 
lengths up to 16 spins. Our results are nicely reproduced by a simple perturbation theory. 
This perturbative approximation is identical to the result of the fermionic effective theory^^ 
when K = TT. We have shown that one needs only to know first-order perturbation theory 
to derive it in a satisfactory manner. For arbitrary field orientation and arbitrary chain 
length we have observed that the perturbative behaviour holds. We thus infer that the 
thermodynamic limit will be also described by the very same approximation. Our ab initio 
results confirm some aspects of the effective theories that have been applied to the spin-1 
chain. 
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The field spfitting of the Haldane gaps in figs. 2,3,4 is that found in experiments 
on NENP^~^^. Both ESR and neutron scattering experiments have observed the same 
sphtting of the Haldane gap. Below the critical field there is clearly a very good agree- 
ment between theory and experiments. At large field the asymptotic behaviour of the 
wavefunctions leads to a polarization of the modes which is similar to experimental data. 

In the future it would be interesting to obtain the magnetization curves in a realistic 
spin-1 following the lines of ref. 19. Another possibility would be to investigate in detail 
the dynamical properties under a field as has been done recently in the zero-field case^°'^^. 
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FIGURE CAPTIONS: 



Figure 1: 

The coefficient kq as a function of the anisotropy D up to D=0.5. The raw data com- 
ing from the finite chain calculation are plotted as open symbols: from bottom to top, 
N=4,6,8,10,12,14,16 and 18. We have performed an extrapolation to the thermodynamic 
limit using Shanks algorithm. The corresponding results are plotted as filled diamonds. 

Figure 2: 

The three gaps in units of J as function of the magnetic field applied along the Z axis. 
The points are results from Lanczos for a 16 spins chain. The solid line is the result 
of perturbation theory eq. (11-14). The dashed lines are the asymptotes of the hyperbola 
(11-14). Note the crossing between the x and z modes. The two points on the right that 
deviate seriously from the perturbative curves are in fact the energies of a state that do 
not belong to the Haldane triplet and that has crossed the upper members of the triplet. 
For the triplet state the deviations always stay small. 

Figure 3: 

The three gaps in units of J as function of the magnetic field applied along the X axis. 
The symbols have the same meaning as in fig. 2. There is no crossing of levels before the 
critical field. 
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Figure 4: 

The three gaps in units of J as function of the magnetic field apphed along the Y axis. 
The symbols have the same meaning as in fig. 2. Note the crossing between the x and y 
modes. 

Figure 5: 

The three gaps in units of J as function of the magnetic field applied along a direction 
close to the Z axis, defined by polar angles 6 = 5 and </> = 45 (in degrees). There is now an 
avoided crossing between z and x modes since there is no longer a discrete symmetry to 
allow a degeneracy. This is generic behaviour since the field is not in a symmetry plane. 

Figure 6: 

The three gaps in units of J as function of the magnetic field applied along a direction 
close to the Y axis {9 = 85 and (j) = 80). The avoided crossing takes place between x and 
y modes. 
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